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We consider a simple direct-sequence model of spread- spectrum 
communications over a bandwidth W using signals of duration T. In 
addition to the dimension n = 2WT of signal space, the other 
parameters of interest are the number (u + 1) of simultaneous users 
of the system, the error rate Pe(u), and the number M of subscribers, 
or potential users. We investigate the relationships between these 
parameters, and, in particular, study the validity of the usual Gaus- 
sian approximation often used to compute Pe (u). Basically, we con- 
clude that, ifM is larger than n 2 /2, then the Gaussian approximation 
is not a guaranteed error rate for (u + 1) users, but rather an average 
over all possible (u + 1) users. IfM is somewhat less than this number 
(the exact value is not known), codes can be assigned so that a 
uniform performance at least as good as Pe(u) can be obtained, 
where, again, Pe(u) is calculated from the Gaussian approximation. 
Uniform performance guarantees are given for any value of M, but 
they (for M large) permit fewer simultaneous users that the Gaussian 
approximation predicts. These bounds explicitly use the maximum 
cross- correlation between the signals of the different subscribers. 
This quantity played no role in the Gaussian approximation. 

I. INTRODUCTION 

In general, spread-spectrum communications refers to a class of 
modulation methods by which the information-bearing signal is trans- 
mitted via a modulated signal having much greater bandwidth. Two 
common methods are used to accomplish the spreading. In one method, 
direct-sequence modulation, the information signal is multiplied by a 
rapidly varying waveform. This waveform, which the receiver is re- 
quired to know, may be thought of as having a pseudo-random char- 
acter. For practical reasons, it has finite duration and is repeated in 
time. We refer to any particular suitable waveform, or a collection of 
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such, as a code. The other method for spectrum spreading is frequency 
hopping. In this case, the available transmission band is divided into 
a large number of disjoint frequency intervals, and the information is 
conveyed by hopping from one such frequency to another. The infor- 
mation may be transmitted by phase-shift-keying each frequency or 
by using a particular set of frequencies for a particular symbol, etc. 
Again the code for frequency hopping must be known to the receiver. 

The initial motivation for introducing such a scheme appears to be 
in its military use as an anti-jamming device. The jammer, not knowing 
the transmitter's code for spectrum spreading, must thus blanket all 
codes. Most of the jammer's power is wasted in codes that are orthog- 
onal to the one in actual use. 

An additional application of more commercial interest was intro- 
duced by Costas. 1 His idea was to use spread spectrum as a way to 
make a large bandwidth, W, available as a communication resource to 
many potential users without preassigning frequency divided channels 
(fdm) (and thus overlimiting the number of potential users) and 
without having a dynamic assignment of fdm (thus incurring the need 
and cost of external control). In modern work, this is usually accom- 
plished, or imagined to be accomplished, by assigning "almost orthog- 
onal" codes, or code vectors, to different users as a means to limit the 
mutual interference between users. 

Very recently, attention has been drawn to spread spectrum as a 
possible modulation method for cellular mobile radio systems. 2 Our 
interest was drawn to this area by Henry's subsequent criticism 3 of the 
analysis of Cooper and Nettleton. 2 While the particular question in 
this controversy appears to have been resolved (in Henry's favor), our 
own survey of the situation has brought out some deeper questions 
relating to assumptions made in the analysis by Cooper, Nettleton, 
and Henry, and often made elsewhere. Specifically, Up service is often 
paid to making the codes approximately orthogonal. Yet when the 
performance analysis is finally made for these digital systems, one 
typical user is considered and all other users which are simultaneously 
using the channel are treated as interfering Gaussian noise having 
uniform power spectrum over the band of interest. Nowhere does any 
measure of the approximate orthogonality of the signals of different 
users enter the performance estimates. Our objective, then, is to 
examine what validity can be given to performance curves calculated 
using the Gaussian approximation and what relation, if any, this 
approximation has to the idea of approximate orthogonality of different 
users' signals. 

To gain some insight into these questions, we examine a simple 
direct-sequence system designed to transmit binary data and give 
upper bounds or the error rate Pe(u) that a user will experience when 
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there are u other users on the channel. Assuming that the bounds yield 
a good description, we can reach some easily stated conclusions. 

Let there be M subscribers or potential users of the system and let 
there be (u + 1) simultaneous users, u + 1 < M. Also denote the 
calculated error rate by Pe(w). Then, by a random coding argument, 
we conclude for M = u + 1 that there is a code such that the (u + 1) 
users each have error rate Pe(w), where Pe(w) is calculated via the 
Gaussian approximation. This, however, does not take into account 
the possibility that the u + 1 users may be selected from M > u + 1 
subscribers. The random coding argument also suggests that, in fact, 
no limit need be placed on M if Pe(u) is estimated from the Gaussian 
approximation. However, the interpretation given to Pe(«) immedi- 
ately changes. It is no longer an error rate for each of (u + 1) users, 
but is an average error rate where the average is taken over all ways 
(u + 1) users are selected from the subscriber population. 

Next, the question of a guaranteed error rate for any (u + 1) users 
is taken up. That is, we present a performance bound valid for each of 
the {u + 1) users which is independent of how the (u + 1) users are 
selected from the population of size M. The results show that if, in our 
n-dimensional signal space, we are packing so many unit energy signal 
vectors (corresponding to different subscribers) that the cross-corre- 
lations (cosines between vectors) are required to be as large as 
(t/n) y2 in magnitude, where t > 1 is simply a convenient parameter, 
then the number of simultaneous users is reduced by a factor of t 
compared to what the Gaussian approximation would predict. Finally, 
recent bounds by Kabatyovskii and Levenshtein for sphere-packing 
problems are applied to give upper bounds on M so that the (t/n ) 1/2 
bound on the cross-correlations can be met. For t= 1, the upper bound 
is n 2 /2 vectors, which, for n large, is quite generous. How closely this 
can be approached is not known. 

II. MODEL AND ANALYSIS 

We consider the following simple model of direct-sequence, binary, 
spread-spectrum communication. We have a bandwidth W over which 
u + 1 independent users simultaneously communicate binary infor- 
mation with a central station at the rate R = 1/T b/s. The individual 
information rates are small compared with the available bandwidth; 
thus, TW » 1. There are M subscribers or potential users of the 
system and M may be large compared with u + 1, the maximum 
number of simultaneous users allowed. Each potential user is perma- 
nently assigned a coded carrier (or code vector). The ith subscriber's 
carrier Ci(t) is written as 

d(t)- £ xji%(t), a) 

A-l 
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where fait) are orthonormal basis functions in the n = 2WT dimen- 
sional space of functions approximately limited to W Hz in bandwidth 
and to T s in duration. 

The entire system operates synchronously and every T seconds the 
ith user square-wave modulates his carrier by ±1, the independent, 
identically distributed (i • i • d) binary data that he wishes to transmit. 

We may consider the n real number x% ] , k = 1, • • • , n forming a 
vector x (I) in real n space and the collection of vectors {x (,) } f belonging 
to the different subscribers is sometimes called a code. We only 
consider equal energy codes and set 

f T 
E = cHt) dt= % (*i°) • (2) 

Jo 



k-i 



In the analysis, we distinguish the user whose performance we will 
be interested in by the subscript i = 1; the other users are designated 
by i = 2, • • •, u + 1. Thus, in a typical T-second interval the received 
signal will be a time translate of 

u+l 

b iCl (t) + J biCi(t), (3) 

i=2 

the bi being the binary data, i«i-d for each user and also between 
users. f We assume correlation detection of (3), and thus the receiver 
bases its decision of b\ on the sign of 

J" T u+i r T u+i 

c\(t) dt + £ ci(t)a(t) dt = Ebi + E % b iPu . (4) 
''" 2 Jo '- 2 

In (4) we have introduced the normalized cross-correlation, 

•T 

d(t)cj(t) dt 



I 



Pij=- = 2 xPxP. (5) 

^ A-l 

This equals the cosine of the angle between x (,) and x^' in ra-space. 
Thus, when | py | is small, the vectors x, and x, are almost orthogonal. 
Assume that b\ = 1 in (4). Then the probability of error, Pe, is 



Pe = Pr 



u+l I 

2 btpu<0 . 



1+ Y bi P u<0\. (6) 



f Such a system where the "spreading" function is modulated by the data is usually 
termed a direct-sequence system. Frequency hopping is another spread-spectrum tech- 
nique. The relative practicality of the two techniques depends on particular circum- 
stances. 
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This depends on the distribution of the noise-like quantity 

q = S b iPu , (7) 

i-2 

where (assuming unit energy) 

P„= l xVxfP. (8) 

A-l 

First let us consider not a specific code, but rather an average over all 
possible codes assuming all the unit vectors x (,) are uniformly and 
independently distributed over the unit sphere.f In this case, note 
(letting < ) denote an average) 

</>h)-- i+h (9) 

n 

which follows from noting that, for a unit vector, the sum of the 
squares of its n direction cosines add to 1, while each must have the 
same average. 

Using the random coding assumption, the moment-generating func- 
tion M Q (s) of q may be shown to be 



M g (s) = <e 89 > 

- (e sv ) u = 



f e av (l - v 2 ) n ~ 3/2 dv\ , (10) 



B(l/2, (n - l)/2) J_j 
where, in (10), B(l/2, (n - l)/2) is the beta function and 

p / V ) 1 (1 _ „2 )( n-3)/2 (11) 

PK ' B(l/2, (n-l)/2) 

is the probability density of any direction cosine v of a vector uniformly 
distributed over a sphere in n-dimensions. A saddle point evaluation 
of (10) yields, for n » 3 and s/n small, 

M v (s) « e s2/2n . (12) 

Since the Chernoff bound 4 states that, for any random variable q, 

Pr[q<A]<e-* A M g {s) for any s < 0, (13) 

then 

Pe = Pr[q < -1] < e 8 M q (s) « e {il,2u/2n)+s) . (14) 



f We could also assume that act" = ±1, i -i -d, and obtain similar numbers. 
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Optimizing the inequality in (12) over s yields + s opt = —n/u, and thus 

Pe(u) < e- (n/2u) , (15) 

where Pe(«) is the error rate for (u + 1) simultaneous users. 

We begin our discussion of (15) by reconsidering the performance 
question from a different point of view. Assume that the interfering 
power resulting from users 2 through u + 1 is uniformly distributed 
over the band W. If each user has power P, then the one-sided power 
spectral density No thus obtained is 

uP uE 
N„ = - = — . (16) 

Further assume this noise is Gaussian. As is well known, the error rate 
for antipodal signals, each of energy E (this is what our transmitters 
are using), is, in white noise of spectral-density No, given by 

(17) 
9<*)--L f e~ y2/2 dy. 



The bound in (17) is exponentially correct. Using (16), the Gaussian 
approximation yields 

Pe(w) < e - {TW/u) = e- {n/2u) , (18) 

where we have introduced the dimension n = 2WT. This is precisely 
the same as the random coding bound (15). We use the random coding 
argument to interpret the result of the standard Gaussian approxima- 
tion. The first interpretation is that there must be a code of (u + 1) 
vectors so that each of the (u + 1) users has error rate Pe(u).* A second 
interpretation is that, to achieve Pe(w), no limit need be placed on the 
number M of subscribers. They can, in fact, be assigned codes at 
random. The average error rate a user sees with u other users present 
is (15). But (15) clearly then refers to a Pe(« ) averaged over all possible 
combinations of u users; it is not an error rate that can be met for any 
set of u other users. Some combinations of u + 1 users will give very 
bad error rates. 

To see what the guaranteed level of performance can be for any 
u + 1 users selected out of a subscriber population of size M, reconsider 



* The fact that s opt = — n/u justifies treating s/n small in (10). 

* More precisely, one can guarantee that a code exists so that a fraction (1 - (1/a)) 
of the users has an error rate no larger than a-Pe(u). 
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(6) with the pu fixed. Then the standard Chernoft bound yields 



Pe < exp 



«+ 1 v-li 



-1(1* 



(19) 



Suppose we have designed our code so that, out of the entire population 
of M subscribers, 

maX |p,y| = Pmax- 

Then from (19), 

Pe<e ,_1/2 " p ™.'. (20) 

If Pmax = l/n, the Gaussian approximation again results. In random 
coding (pi) =l/n, but nothing was said about p max and thus nothing 
could be said about a guaranteed error rate. 

Given M unit vectors in n space, how small can p max be? One result 
in this direction is due to Welch. 3 He states that 



Pmax ^ 



n 



(21) 



M-l 
If n is large and M large compared to n, then (21) already states 

2 1 t 

Pmax > - • 

n 

Rephrasing our question, given the dimension n and p ma x, how large 
can M be? That is, how many vectors u, can we put on the unit sphere 
so that 

Ia/I-Iuj-UsI^jd.*.? (22) 

This sphere-packing problem is different from the conventional one 
which requires 

pij = u,-u 7 <p ma x. (23) 

The number of vectors M will be much smaller under condition (22) 
than under condition (23). Luckily, our sphere-packing problem (22) 
was one of the packing problems recently considered by Kabatyonskii 
and Levenshtein. 6 Precise values of M are not known, but upper 
bounds are. 

An extension of their work provides us with the following. For 

A-0,1,2, ..., 



* If (n + 1) vectors are the vertices of the regular simplex in n-space, then pL. = 
l/n 2 . Then (21) is exact for M = n and M = n + 1. 
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(1-pLx) . , .. (n/2) k 

M * (2k + l)-(n + 2k) P lJ n + 2k) OM k > m 

provided the denominator in (24) is positive. + In (24), we have used 
the notation 

(a) k = a(a + 1) • • • (a + k - 1) 

(a)o = 1. (25) 

It is useful to write piLx as 

Pma* = ~ • (26) 

n 

If we set t = k = integer and use the same k in (24), and further assume 
k 2 <K n, we have approximately 

1 n k+1 

Mrs . (27) 

(A + 1) (2k - 1)!! l ' 

Equation (24) is plotted in Fig. 1 for n = 100 as a function of t. The 
part of the curve for smallest t uses k = 0, the next set of t values uses 
k = 1, and so on. 

In general, the value of the upper bound at t = 1 (where the Gaussian 
assumption agrees with random coding) is given by n 2 /2. How closely 
this upper bound can be achieved is not clear. If, however, t = 4, a 
result of van Lint and interpreted for the present situation by Welch, 5 
indicates that at least 10 4 binary (±1) waveforms are available, com- 
pared with the bound of 4 X 10 5 . The point is indicated by the small 
circle in Fig. 1. 

III. CONCLUSIONS 

Consideration of our simplified model has provided the following 
insights. Treating many other users as background Gaussian noise for 
a particular user is a good approximation in that a code can be found 
which (at least approximately) provides the calculated error rate for 
each simultaneous user. This is no bargain, however, for the orthogonal 
code would be even better in performance in the present problem, 
permitting more simultaneous users.* Evaluated as a pure modulation 



f The bound (24) is valid for a real vector space. For a complex space, the following 
larger bound applies (p = p m *x): 

•Mcomplex^- 



k + 1 - (n + k)p 2 M 



* More explicitly, (18) states that, for a small error rate, we require 2u ■« n = 2 WT, 
or u <sc WT = W/R. However, in this model, fdm with double-sideband modulation 
allows about W/R users with zero error rate. 

2020 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1979 




Fig. 1 — Upper bound on number of signals M vs t 
is known to be achievable using binary codes. 



npL for n = 100. The point A 



method for accommodating (u + 1) known and fixed users in a white- 
noise channel, spread spectrum has little to offer over ssb-fdm. The 
only advantage for spread spectrum on the type of channel that we 
have considered would be as a multi-access scheme for M potential 
users. + Selecting codes at random for this situation still permits us to 
use the Gaussian approximation with u "other" simultaneous users to 
calculate an error rate Pe(w), but only in an average sense. Namely, it 
is an average over all possible (u + 1) users selected out of the M 
subscribers. To guarantee a level of performance for any (u + 1) users, 
the departure from strict orthogonality must not be too severe, and 
this puts a definite restriction on the number M of potential users. 



+ The potential absence of channel assignment for spread spectrum is the basis of the 
Cooper-Nettleton (Ref. 2) proposal for spread spectrum for mobile radio. Other consid- 
erations may make spread spectrum an attractive alternative. For example, it can 
provide frequency diversity for a frequency-hopped system when there is fading. Such 
a scheme has been proposed by Goodman et al. (Ref. 7) for mobile radio, modifying 
Viterbi's (Ref. 8) proposal for satellites. Also in satellite systems, the Doppler shift can 
be large compared to data rates for individual users. Viterbi (Ref. 8) has suggested that 
spread spectrum would not need the large guard bands between channels that fdm 
would require. 
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While upper bounds on M were given in the text, the exact value is not 
known, nor, in general, did we discuss explicit construction of subop- 
timum codes. 
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APPENDIX 
Derivation of Eq. (24) 

A brief outline is presented to guide the reader who wishes to 
rederive (24) from the results of Ref. 6. 

We first summarize the relevant results of Ref. 6, given on pages 10 
and 11 of that work. We will be concerned with the interval — 1 < t < 
1 and expansions in Jacobi polynomials Pfit), i denoting the degree 
of the polynomial. The parameters a and /? depend on the dimensional- 
ity n and type (real or complex) of space that we are considering. For 
a real space a = (n — 3)/2, /? = — 1/2 (n > 3), while for a complex 
space a = n — 2, (S = (n > 2). Let s be a real number — 1 < s < 1 and 
denote by R (a, /?, s) the set of polynomials 

fit) = 2 fiPfit) (28) 

i-O 

of degree / = 1, 2, • • • such that 

(i) /i>0, i = 0,...,k where f > 0. (29a) 

(ii) f(t)£Q for -1<*<s. (29b) 

Then the maximum number M of unit vectors u, in /i-space such that 

|u,-u,|<p 

satisfies 

Mrs inf ^. (30) 

The evaluation of f(l)/fo for any particular allowed f(t) further upper 
bounds (30). We choose the polynomials 

f(t) = (t-s)(t+l) k k = 0, 1, .... (31) 

The verification of (29b) is thus trivial, as is the evaluation of /(l). The 
verification of (29a) as well as the evaluation of fo is a direct calculation. 
We evaluate /j using the orthogonality properties of the Jacobi poly- 
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nomials Pf{t) with respect to the weight function w(t) = (1 - t) a (l 
+ t) fi . The required normalization may be found in Ref. 9, p. 262, 
formula 1. If we rewrite (31) as 

f{t) = (t+ l) k+1 - (1 + s)(t + l) k , (32) 

then the integrals needed to evaluate f may be calculated from Ref. 8, 
p. 263, formula 3. One then directly verifies that f, > {i > 0) is positive 
whenever f is, and the indicated evaluation of f results in (24) for the 
real case. 
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